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The existence of a normal state spectral gap in underdoped cuprates raises important questions about the 
associated superconducting phase. For example, how does this pseudogap evolve into its below Tc counterpart? 
In this paper we characterize this unusual superconductor by investigating the nature of the "residual" pseudogap 
below Tc and, find that it leads to an important distinction between the superconducting excitation gap and order 
parameter Our approach is based on a conserving diagrammatic BCS Bose-Einstein crossover theory which 
yields the precise BCS result in weak coupling at any T < Tc and reproduces Leggett's results in the T = 
limit. We explore the resulting experimental implications. 
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Pseudogap properties, associated with the unusual normal 
state of the underdoped high temperature superconductors, 
have received considerable attention in the literature. From an 
experimental perspective the relationship (if any) between the 
pseudo- and superconducting gaps has not beeruunambigu- 
ously clarified. Angle-ijesplved photoemissionBcl (ARPES), 
and other measurementgja on the underdoped cuprates indi- 
cate that the normal state excitation or, equivalently, pseudo- 
gap above Tc evolves smoothly into the excitation gap at and 
below Tc- It is unlikely that a fully developed pseudogap will 
abruptly disappear as the temperature falls below Tc, but pre- 
cisely how it connects with the superconducting order param- 
eter is not obvious. While there are scenarios in which the 
pseudogap and superconducting gaps are believed to be inter- 
relatedjo here we point out a quantitative relationship for one 
widely discussed scenario and characterize key features of the 
resulting unusual superconducting state. In this way we sug- 
gest experimental tests which may distinguish one pseudogap 
model from another. 

A large class of pseudogap scenarios for the underdoped 
cuprates associate, this phase with some form of precursor 
superconductivity!], often in the context of a (normal) state 
intermediate between that of the free fermioM of the BCS 
and bound pairs of the Bose-Einstein regimesB. These BCS 
Bose-Einstdn crossover theories were originally formulated 
by LeggetJj to address the nature of the superconducting 
ground state. j-There has been considerable attention paid to 
this approachll2rt3, primarily at and above Tc- Our goal here 
is to establish a crossover theory in the regime < T < Tc, 
which is consistent with three important criteria. These in- 
volve simultaneously (i) satisfying the law of particle conser- 
vation, (ii) establishing consistency with the precise BCS re- 
sult in weak coupling, and (iii) establishing consistency with 
the formulationpoJjJJef. ^ for the ground state. Of jthese three 
criteria, the firsttaEJ and second, as well as thirdll^M have not 
necessarily been satisfied in earlier work. In this process we 
determine the counterpart to the pseudogap below Tc, and its 
experimental implications. 

In this paper we build on previous work0~lli, based on a 
particular diagrammatic version of these crossover theories. 
For definiteness, we take a simple model of 3D fermions 
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which interact via a short range, separable pairing interac- 
tion with s-wave symmetry V^,,/ = —\g\^Pkfk'- It should 
be stressed that we have previously demonstrated that our re- 
sults (above Tc) for this isotropic model remain qualitatively 
similar when applied tOpa,quasi-2D lattice model with attrac- 
tive c?-wave interactionsllj. While, the latter is more suitable 
for describing the superconducting state of the cuprates, the 
general physics we discuss here is presented more clearly, 
without the complexity of d-wave pairing. Our diagrammatic 
scheme is ba^ed on the "pairing approximatiorClpf Kadanoff 
and MartincJ, subsequently extended by Pattoncil, which will 
be shown below to satisfy the three criteria discussed above. 
Following these references, one arrives at the following com- 
plete set of equations 

nK) = G-HK)~G-\K) 

^J2*{Q)Go{Q-K)^l_^^,, (la) 
Q 

5= [l + .9X(Q)]t(Q) , (lb) 

x(Q)=5]G(X)Go(Q-if)^Lq/2, (Ic) 

K 

n = 2^G(if), (Id) 
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which self-consistently determine both the Green's function 
G{K) and the T-matrix t{Q). Equation ( p^ ) is associated with 
particle conservation, criterion (i). For brevity, in Eqs. (1) we 
have used a four-momentum notation K = {k;iu}) and Q = 
{q; i^l), where uj/fl are odd/even Matsubara frequencies. The 
bare Green's function is given by Go{K) — {iuj ~ ■^k)^^, 
with = k"^ /2m — fi and n is the particle number density. 
Here S(i^) is the self-energy and x{Q) the pair susceptibil- 
ity. In the weak coupling limit, Eqs. (1) can be regarded as a 
T-matrix formulation of the generalized BCS theory with 

tsciQ) - --^KQ) , (2) 

where Age is the superconducting order parameter, and the 
Dirac-delta function guarantees the factorization of the two- 
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particle correlation function in a manner consistent with off- 
diagonal long range order From Eq. (|l|), the corresponding 
BCS self-energy is given by 
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With (||) and (||), Eqs. (1) yield the usual BCS gap equation for 
Asc- As the coupling strength g is increased, the role of pair 
fluctuations (representing the mean square deviation of the 
pairing field from its average value |Asc|) becomes increas- 
ingly important and additional contributions to the T-matrix 
need to be appended to Eqs. (^. These effects are precisely 
those needed to describe pseudogap phenomena above Tc- 
We write the T-matrix below Tr as 



tiQ) = tsc{Q)+tpg{Q) , 



(4) 



where the "singular" tsc, given by (g), accounts for the con- 
densate of Cooper pairs, while the "regular" tpg describes pair 
fluctuations associated with the pseudogap. Inserting into 
Eq. (p^, along with and using the filtering property of the 
Dirac-delta function, one obtains 



tpg {Q) 
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and 



l+.gx(0) = 
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This last equation is the self-consistent gap equation. More- 
over, we see from the above two equations that the pseudogap 
component of the T-matrix, tpg{Q), isJnighly peaked about 
the origin, with a divergence at Q = Ocj. The self-energy of 
Eq. ( [la[ ) may be decomposed into two contributions 
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In evaluating the pseudogap contribution to the self-energy, 
detailed numerical calculationstS show that the main contri- 
bution to the Q sum comes from the the small Q region so 
that J^pgiK) = Go{-K) J2q tpg{Q)+5^{K), where JS, 
in the momentum and frequency range of interest, is much 
smaller than the leading BCS-like part, and can be ignored in 
what foUowscj. Thus, E^g can be well approximated by the 
BCS-like form 
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where the pseudogap amplitude within the superconducting 



state, Apg, is defined as 



Q 



tpgiQ) 



(9) 



Note that, although |Asc| satisfies an equation similar to (^), 
there is an important distinction between the two energy gaps. 
Asc is a complex order parameter which represents the mean 



value of the pairing field. By contrast, the pseudogap param- 
eter Apg is a positive definite quantity which describes the 
(incoherent) fluctuations of the pairing field about its mean 
value. 

It foUows from Eqs. (0)-(|l) that 
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The central results for the superconducting gap equa- 
tion, number density and pseudogap below Tc follow from 
Eqs. ( 1 ,p|Jq,p|,[10|) and can be summarized as 
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where Ek = (Cfe + A^(/3^)^/^. Note that the total excita- 
tion gap A and the chemical potential ca n be obta ined from 
first two Eqs. dill). Moreover, while Eqs. (Ha lit) coincide 
formally with the corresponding (weak coupling) BCS equa- 
tions, here these equations are valid for arbitrary coupling and 
any T < Tc. It should be stressed that Equation (11c), which 
determines the precise decomposition of A into Age and Apg, 
is crucial and contains much of the central new physics of this 
paper 

The simplicity of these equations derives directly froipJhe. 
diagrammatic scheme of Eqs. (1); alternative schemesll3'E£l 
will not produce this standard form, nor will they lead to the 
BCS limit in the weak coupling case. Moreover, as will be 
seen below, w hen T 0, the pseudogap Apg vanishes ani 



Eqs. ( |1 lalJl lb| ) coincide precisely with those used by Leggeti 
in his T = BCS Bose-Einstein crossover theory. Finally, as 
T ^ Tc from below, the equations satisfied by Tc, Apg andjj. 
can be seen to be identical to their counterparts found earlierlla 
when Tc is approached from the normal state. 

Physically, the pseudogap below Tc can be interpreted as an 
extra contribution to the excitation gap, reflecting the fact that 
at moderate and large g, additional energy is needed to over- 
come the residual attraction between excited fermion pairs in 
order to produce fermionic-like Bogoliubov quasi-particles. 
In the bosonic limit, it becomes progressively more difficult 
to break up these pairs and the energy Apg increases accord- 
ingly. 

In Fig. [l|a are plotted Apg, Age and A, as a function of tem- 
perature, obtained from a numerical solution of Eqs. ([TTI). We 
choose for illustrative purposes three representative values for 
g/gc — 0.7, 0.85 and 1.0, (all of which lead to positive chem- 
ical potential), corresponding, respectively, to a small, inter- 
mediate and large pseudogap parameter at Tc- Here, for defi- 
niteness, we follow Ref. [l^ and take ipk — {I + fc^/fco)^^/^, 
with kg = Akp and define gc ~ —Air/mko to represent the 
critical coupling necessary to form a bound fermion pair in 



2 



vacuum. As can be seen from the figure, with decreasing tem- 
perature, Apg{T) decreases mono tonically, from its maximum 
value at Tc until it essentially vanishesEa at T = 0, while 
Asc(r) and A(r) both increase monotonically. When ap- 
proached from slightly above Tc, there will be a slope discon- 
tinuity in A at Tc, reflecting the related discontinuity in the 
order parameter Age', moreover, as demonstrated in Fig. |l]a, 
at the higher value of g this total gap is almost temperature 
independent. 

The pseudogap is an important measure of the distinction 
between the order parameter, A^^, and the excitation gap A. 
The latter is the quantity deduced in ARPES measurements. 
The former must be directly related to the superfluid density. 
Us, which is strictly zero at and above Tc, but which presum- 
ably also depends in some way on A as well. Moreover, by 
studying rig, (which can be obtained via the London penetra- 
tion depth), it will be made clear that, in principle, all three 
gap parameters can be distinguished experimentally. The su- 
perfluid density may be expressed in tacms of the local (static) 
electromagnetic response kernel K{O jEB 

in in 
ng = ^K{0) = n-—^Pc,c.{0), (12) 

with the current-current correlation function given by 

Po^piQ) = -2e2^Aa(if,i^ + Q)G(i^ + g) 

K 

xAp{K + Q,K)G{K) . (13) 

Here the bare vertex X{K, K + Q) = ^{k + q/2), while the 
renormalized vertex A must be deduced in a manner consis- 
tent with the generalized Ward identity, aaplied here for the 
uniform static case: Q — (g,0), q OEs. It is convenient 
to write A = A + ^A. 
tribution ^A 
identity 



pg ^ 5\sc, where the pseudogap con- 
to the vertex correction follows from the Ward 
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The particle density n, given by Eq. ( |ld[ ), after partial inte- 
gration can be rewritten as n = —(2/3) fe • dG{K) / dk. 
Then, as a result of Dyson's equation, one arrives at the fol- 
lowing general expression 
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Now, inserting Eqs. ( p3| ) and ( p3| ) into Eq. (|l2|) one can see 
that the pseudogap contribution to Us drops out by virtue of 
Eq. (14); we find 



lY.G\K)k.{5A. 



dk 



(16) 



We emphasize that the cancellation of this pseudogap contri- 
bution to the Meissner effect is solely the result of local charge 
conservation. 




"(a) T/To T/Tc (b) 

FIG. 1. (a) Temperature dependence of the excitation gap A, su- 
perconducting gap I Asc| and pseudogap Apg for coupling strengths 
g/gc = 0.7, 0.85 and 1.0. (b) Temperature dependence of the super- 
fluid densities tIj, and rff for the same coupling strengths as in (a). 
The shaded regions emphasize pseudogap effects. 

Following the standard prescription for constructing the 
proper vertex cep-ection corresponding to the superconduct- 
ing self-energyEa one obtains 

SA.ciK + Q,K)^ Al^^lGoi-K - Q)Go{K)\{K + Q) . 

(17) 

Inserting Eqs. (||,0,|7|) into Eq. ([l^, after calculating the Mat- 
subara sum, one arrives at 
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We may write the superfluid density as Ug = A^^ F{A) 
[where the form of the function F can be obtained from 
Eq. (p^]. The same quantity corresponding to a BCS su- 
perconductor with effective gap parameter A is given by 
nf = A^F{A), so that, Us/nf = {Asc/Af < 1. In 
Fig. |l]b we plot the temperature dependence of the normalized 
superfluid density ns/n (solid line) calculated from Eq. (|l8|) 
for the same three representative values g/gc as above. These 
curves are compared (dashed line) with the quantity nf^/n, 
which is a (BCS-like) function only of the excitation gap. For 
sufficiently weak coupling [g/gc ^ 0.7) the two curves are 
indistinguishable. With increasing g the separation between 
the two curves become evident, particularly in the vicinity of 
Tc, whereas at zero temperature there is no difference since 
lis ~ n, independent of the coupling. This comparison thus 
demonstrates how different are these "pseudogap" supercon- 
ductors. The superfluid density reflects most directly the tem- 
perature dependence of Asc, not the excitation gap. 

The existence of residual pairing correlations below Tc will 
affect thermodynamic properties as well. Indee^upon analy- 
sis of data in underdoped cuprates, Loram et altll conjectured 
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that the measured excitation gap squared can be expressed 
as the sum of the squares of a pseudogap and superconduct- 
ing order parameter This pmely^ienomenological analysis 
leads to a similar decompositionc3 of the excitation gap, as 
in Eq. ( [lO| ) However, in contrast to the present work, these 
authors presumed that Apg is temperature independent below 
Tc. 

In summary, in this paper we have demonstrated that, if a 
pseudogap state arises from pairing correlations (fluctuations) 
above Tc, then these pairing fluctuations necessarily persist 
below Tc. These pseudogap systems are unconventional su- 
perconductors, in which pair fluctuations are present all the 
way down to the lowest temperatures. At T = these fluc- 
tuations (or Apg) vanish. A key manifestation of the "su- 
perconducting pseudogap" is in the nature of the excitation 
gap (A), which differs significantly from the superconducting 
order parameter Age, as A^ = A^^, + A^^. At a physical 
level we view Apg as reflecting an additional energy associ- 
ated with the attractive interaction, which must be overcome 
in order to create fermionic-like Bogoliubov quasi-particles. 
In this way, the excitations from the condensate in a BCS 
Bose-Einstein crossover theory can be viewed as intermediate 
between the (free) fermionic Bogoliubov quasi-particles of the 
BCS limit and the (bound) bosonic pairs in the Bose-Einstein 
regime. It should be jStressed that our previous work on d- 
wave superconductorslij reinforces the claim that the physics 
presented here for the s-wave case, is not qualitatively sensi- 
tive to the symmetry of the pairing interaction. 

Experimentally, verification of this pseudogap scenario (for 
the underdoped cuprate superconductors) involves establish- 
ing the relation between A and Age- Measurements of Ug 
and A separately are possible (through penetration depth and 
ARPES experiments). Even more promising may be tunnel- 
ing spectroscopy measurements of high Tc superconductor- 
insulator-superconductor junctions in which the Josephson 
and quasiparticle current data can be simultaneously used to 
extract A and Age- 
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